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Abstract. Reconstruction of jets in high-energy heavy-ion collisions is challenging due to the 
'^r large and fluctuating background coming from the underlying event. We report results on full 

Cn jet reconstruction, obtained from data collected in 2011 by the ALICE detector at LHC for 

^_^ Pb-Pb collisions at -^snn = 2.76 TeV. The analysis makes use of the tracking system and 

K> the electromagnetic calorimeter. Signal jets, which come from hard scattered partons, are 

QJ reconstructed using the anti-fcT jet finder algorithm. The average background is subtracted on 

' , a jet-by-jet basis to reduce the contribution to the jet reconstructed energy coming from the 

O underlying event. The jet spectrum is corrected to account for fluctuations in the background 

^ momentum density and detector effects through unfolding. 
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OO 1. Introduction 

/^ QCD jets are produced in high-energy particle cohisions as a result of the fragmentation of 

\^ a high momentum scattered parton. Experimentally they are reconstructed using a well- 

.• defined algorithm, which acts as a working definition of a jet, to be used consistently also 

^^ in phenomenological models. In heavy-ion collisions, hard scattered partons are produced in 

C^ the early stages of the collision, so that they propagate through (and potentially are affected 

^\~l by) the hot and dense nuclear medium. The interactions suffered by the parton can result in 

t> energy loss, widening and/or complete absorption of jets. These phenomena go under the name 

k> of "jet quenching". RHIC and LHC experiments have already collected convincing evidence of 

Vh jet quenching in a number of measurements, such as the suppression of high-px particles relative 
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to a pp baseline [H El 13 HI , hadron-hadron correlations [5", "6] and jet-jet correlations at high- 
Pt [3 IS]- However, performing a full jet reconstruction in the low and intermediate px region 
is particularly challenging due to the overwhelming soft particle background coming from the 
underlying event (UE). ALICE has recently reported on a first detailed study of the background 
for jet reconstruction in the heavy-ion environment ^ . 

In these proceedings the analysis techniques utilized to fully reconstruct jets in Pb-Pb 
collisions with the ALICE experiment are presented. The results shown here come from data 
collected by ALICE in fall 2011 at ^/si^ = 2.76 TeV. 

2. Inputs to the jet finder 

The ALICE tracking system benefits of the Inner Tracking System, a six-layer silicon detector 
which provides a precise measurement of the primary vertex together with the first points of the 
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(a) R = 0.2. 
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(b) R = 0.3. 



Figure 1: Raw jet pT spectra in Pb-Pb collisions at ^snn = 2.76 TeV in the 0-10% centrality 
range at mid-rapidity: inclusive is shown in full circles and with increasing minimum leading 
hadron p^ requirement in the other symbols (color online). Jets are reconstructed using the 
anti-A:T algorithm with two different resolution parameters, R = 0.2 (left) and R = 0.3 (right). 



tracks, and of a large Time Projection Chamber [lOj. Tracks are reconstructed at mid-rapidity 
{\r]\ < 0.9) and in full azimuth. The ALICE electromagnetic calorimeter [TT] is a Pb-scintillator 
sampling calorimeter, which covers mid-rapidity (|r/| < 0.7) and partial azimuth (Ac/? = 100°). 
The EMCal measures photons, e.g. from vr'^ decays, which are included in the jet finder input. 
The shift of the jet energy scale due to unreconstructed particles, such as K^ and neutrons, has 
been studied in pp simulations and accounted for in the final result. A full description of the 
ALICE experiment is available in Ref. |10| . 



3. Jet finding and average background 

The anti-/cT [IS] jet finding algorithm has been employed in its Fast Jet [l3j implementation. 
This sequential recombination algorithm has the advantage of being "soft-resilent" , namely it 
is little affected by the soft background JT^■ Anti-fcT jets are pretty regular, cone-like shaped 
around some high-pT particle. The uncorrelated energy density from soft processes is large in 
heavy-ion events and needs to be subtracted on an event-by-event basis. The main problem in 
the estimation of the average background is the separation of the UE from the hard scattering. 
The method used here has been proposed in Ref. p^. The average background p is calculated, 
event-by-event, as the median of the pT density (jet pT over jet area) of the fex algorithm [16] 



reconstructed jets. A jet-by-jet subtraction is performed: Px let ~ Pt^]"" ~P^ ^jet ; where Pt^?°" 
and Ajet are respectively the transverse momentum and the area of the jet. Combinatorial jets, 
i.e. jets reconstructed out of the soft background, are efficiently removed by requiring a minimum 
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Figure 2: (5pT distributions for 
Pb-Pb collisions at ^snn = 
2.76 TeV in the 0-10% central- 
ity class. Two different meth- 
ods have been implemented: "ran- 
dom cones", shown with full sym- 
bols, and single particle embedding, 
shown with open symbols (see text 
for details); the distributions for 
R = 0.2 are narrower w.r.t. R = 0.3 
(color online). 



leading hadron px- Figure [T] shows the anti-fcx raw jet pT spectra obtained applying different 
minimum leading hadron p^ requirement, from (inclusive) to 10 GeV/c, for two different jet 
resolution parameters R = ^^/Arf + Aip^. 

4. Unfolding 

The spectra shown in Fig. [T] are not directly comparable with model predictions. In fact the 
measured values of the observable, the jet px, are subject to random fluctuations, due to region- 
to-region differences in the background momentum density and to the detector response. This 
means that each observation is characterized by a true (and unknown) value Px^ef ^^^ ^y ^ 
measured value Pxlet- '^^^ histograms for p^^et ^^^ ^xlet ^^^ related by a convolution through 
the response matrix RMtot = RMdet x R-Mbkg, where RMjct parametrizes the detector reponse 
whereas RMbkg describes background fluctuations. Unfolding is the numerical procedure that 
allows to get back the true distribution given the measured one and the response matrix RMtot • 

Background fluctuations have been estimated in two different ways [9] , namely using random 
cones (scalar sum of the pT of all particles found in a cone randomly placed in the event) 
and single particle embedding (with the anti-fcx algorithm). The residual px differences due to 
region-to-region fluctuations are calculated as: (5px = PTjet — PT.probe — P ^ ^jetj where px.probe 
is the pt of the embedded probe (px,probe = for random cones). The 5pT distributions, shown 
in Fig. [2| tell us how much the jet px is smeared due to background fluctuations. 

The detector response to jet reconstruction has been studied with pp simulated events, using 
the PYTHIA6 [T7j generator and the GEANT3 |18| transport code. Jets are reconstructed both 
at generator level and at detector level. The generator-level and detector-level jets are matched 
following a geometrical criterion. 



5. Results 

This measurement makes use of the unfolding iterative method proposed by D'Agostini [19], 
which contains elements of Bayesian statistics. In Bayesian unfolding the number of iterations 
plays the role of the regularization parameter. Based on the evolution of the covariance matrix 
and the converging of the solution itself, a number of four iterations has been chosen as 
default. Indeed after three/four iterations the procedure starts to converge. However, with 
more iterations a characteristic fluctuation pattern does appear along with larger variances and 
(anti-)correlations between far bins (indicating under- regularization) . Systematic uncertainties 
have been estimated taking into account variations of the solution for ±1 iterations (~ 5%). 
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Figure 3: On the left: spectra of reconstructed jets in Pb-Pb collisions at -^snn = 2.76 TeV in 
the 0-10% centrality class with two different niinimum leading hadron pT requirements. On the 
right: ratio of the two spectra (solid circles) and comparison with a PYTHIA pp expectation 
(open circles). The bands represent the systematic uncertainty; statistical uncertainties are 
shown with bars, when they are larger than markers (colors online). 



To better understand the effect of the leading hadron requirement in the jet pT spectrum, 
the analysis was performed for PT,ieading > 5, 10 GeV/c. Figures Is] show the two spectra with 
the two different leading hadron thresholds, and the ratio, compared to a PYTHIA expectation. 



6. Conclusions 

We have reported on the analysis techniques utilized to reconstruct jets in the 10% most central 
Pb-Pb collisions at ^/snn = 2.76 TeV recorded by ALICE in 2011. The raw jet px spectra 
are corrected for the average background and biased requiring a minimum leading hadron p^. 
Corrections for background fluctuations and detector effects are applied via Bayesian unfolding. 
The effect of the leading hadron requirement was studied for two different thresholds, 5 GeV/c 
and 10 GeV/c. The ratio between the two corrected spectra is in reasonable agreement with a 
PYTHIA simulation, which indicates a vacuum-like fragmentation of the jet core. 
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